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DIRECT  SIMULATION  MONTE  CARLO  STUDY  OF  H/Hj  AND  H/H2/CO 
MIXTURES  FOR  DIAMOND  CHEMICAL  VAPOR  DEPOSITION 

I.  INTRODUCTION 

Chemical  vapor  deposition  (CVD)  has  become  the  dominant  technique  for  the 
growth  of  high-quality  diamond  films  possessing  exceptional  mechanical,  optical,  ther¬ 
mal,  and  electrical  properties  [1].  Successful  diamond  growth  has  been  demonstrated 
over  a  wide  range  of  conditions  using  microwave  plasma,  hot  filament,  combustion 
flame,  and  plasma  torch  reactors.  These  various  techniques  can  be  classified  as  ei¬ 
ther  diffusion-dominated  or  convection-dominated  on  the  basis  of  their  characteristic 
Peclet  number  [2],  defined  as  Pe  =  uLfD,  where  D  is  the  diffusion  coefficient,  u 
is  the  flow  velocity,  and  L  is  the  length  scale  for  the  transport  of  species  to  the  di¬ 
amond  surface.  The  microwave  plasma  and  hot  filament  reactors  typically  involve 
very  slow  flow  velocities  and  fall  well  within  the  diffusion  dominated  regime  (Pe<Cl), 
while  the  combustion  flame  and  plasma  torch  reactors,  which  rely  on  the  high-speed 
forced  convection  of  gases  to  the  diamond  surface,  are  clearly  convection-dominated 
(Pe>l). 

The  focus  of  this  paper  is  on  the  diffusion-dominated  processes.  In  this  regime, 
convective  gas  flow  can  be  neglected  and  the  problem  can  be  effectively  reduced  to 
that  of  diffusion  through  a  one-dimensional  boundary  layer.  To  study  this  system,  we 
use  the  direct  simulation  Monte  Carlo  (DSMC)  method  [3].  This  approach  allows  us 
to  accurately  extend  our  calculations  over  a  range  of  conditions  from  those  where  a 
Navier-Stokes  continuum  treatment  is  valid  into  a  regime  where  rarefied  gas  effects  are 
important.  Since  DSMC  is  a  particle  based  method,  the  effects  of  molecular  diffusion, 
thermal  diffusion,  viscosity,  and  thermal  conduction  are  all  implicitly  treated. 

The  primary  aim  of  this  work  is  to  determine  the  effects  of  the  operating  condi- 
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tions  on  the  hydrogen-atom  and  thermal  fluxes  to  the  substrate.  A  number  of  studies 
have  shown  the  critical  role  of  hydrogen  atoms  in  both  creating,  through  hydrogen 
abstraction  reactions,  surface  radical  sites  where  growth  species  can  be  added  [4, 
5],  and  in  “annealing”  the  surface  by  promoting  sp^  bonding  among  carbon  atoms 
newly  added  to  the  diamond  lattice  [4,6].  Since  we  are  not  here  concerned  with  the 
specific  mechanisms  of  carbon  addition,  simplifications  can  be  made  concerning  the 
composition  of  the  gas  phase. 

Standard  diamond  CVD  source  gas  recipes  consist  of  a  small  amount  (<  1  %)  of 
methane  in  hydrogen.  The  relatively  small  concentrations  of  methane,  methyl  radi¬ 
cals,  and  other  carbon  species  will  have  negligible  effects  on  the  thermal  conductivity 
of  the  mixture  and  the  diffusion  of  atomic  and  molecular  hydrogen.  Such  mixtures 
can  be  approximated  in  our  present  simulations  as  pure  hydrogen.  Similar  simplifying 
assumptions  can  be  made  for  source  gases  containing  oxygen.  Bachmann  et  al.  [7] 
have  shown  that  source  gas  recipes  containing  oxygen  can  be  categorized  as  leading 
to  diamond  growth,  non-diamond  growth,  or  no  growth  based  solely  on  their  overall 
C-H-0  stoichiometry.  This  independence  of  the  particular  source  gases  employed  sug¬ 
gests  that  a  quasi-thermodynamic  equilibrium  is  established  before  the  gases  reach 
the  diamond  surface.  Consequently,  the  carbon  and  oxygen  in  the  system  will  react 
to  form  as  much  carbon  monoxide  as  possible,  leaving  the  excess  element  to  form 
either  growth  (CH3,  C,  etc.)  or  etching  (OH,  O)  species.  Since  the  highest  quality 
diamond  growth  occurs  for  mixtures  with  only  a  small  excess  of  carbon  to  oxygen, 
it  is  reasonable  to  model  these  systems,  for  the  purpose  of  this  study,  as  H/H2/CO 
mixtures.  Because  of  our  interest  in  diamond  growth  from  ethanol-water  source  gases 
[8-11],  we  have  chosen  an  H/H2/CO  mixture  with  a  stoichiometry  identical  to  that 
obtained  from  a  1:1  ethanol-to- water  mixture. 

The  DSMC  method  and  the  details  of  the  numerical  simulations,  including  the 
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gas  phase  and  surface  chemistry  models,  are  described  in  §11.  This  is  followed  by  a 
presentation  of  the  numerical  results  for  a  range  of  operating  conditions  in  §III.  In  §IV , 
the  DSMC  results  are  analyzed  and  compared  to  analytic  expressions  describing  the 
behavior  of  the  system  in  the  low-  and  high-pressure  limits.  Details  of  the  derivations 
of  the  analytic  forms  are  given  in  the  Appendix.  Finally,  the  implications  of  our 
findings  for  diamond  CVD  are  discussed  and  summarized  in  §V. 

II.  NUMERICAL  MODEL 

In  this  section  we  give  a  brief  overview  of  the  direct  simulation  Monte  Carlo 
(DSMC)  method.  For  a  more  comprehensive  treatment,  see  Bird  [3].  This  is  followed 
by  a  description  of  the  one-dimensional  computational  domain,  the  gas-phase  collision 
model,  and  the  surface  chemistry  mechanism. 

A.  DSMC  method 

The  degree  of  rarefaction  of  a  gas  flow  is  usually  measured  by  the  Knudsen  num¬ 
ber,  Kn  =  A/L,  where  A  is  the  mean  free  path  of  the  gas  molecules  and  L  is  the 
smallest  characteristic  dimension  of  the  system.  The  traditional  Navier-Stokes  treat¬ 
ment  of  gas  dynamics  breaks  down  at  approximately  Kn  =  0.1,  and  must  be  replaced 
by  an  appropriate  molecular  model.  Although  molecular  dynamics  (MD)  is  valid  in 
this  regime,  its  implementation  would  be  prohibitively  expensive.  An  alternative  ap¬ 
proach  is  the  DSMC  method  developed  by  Bird.  DSMC  is  a  particle  method  in  which 
the  molecular  motion  is  modeled  deterministically,  while  the  binary  collisions  between 
molecules  are  treated  statistically.  Although  DSMC  has  not  been  strictly  proven  to 
be  equivalent  to  solving  the  Boltzmann  equation,  the  technique  is  based  on  the  dilute 
gas  and  molecular  chaos  assumptions  used  in  its  derivation.  The  DSMC  method  has 
been  used  primarily  by  the  aerospace  community  in  the  study  of  rarefied  flows  [3], 
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but  more  recently  has  been  applied  to  the  modeling  of  CVD  processes  [12-15].  The 
applicability  of  the  method  is  not  limited  to  rarefied  flows,  but  because  the  computa¬ 
tional  expense  of  the  method  increases  proportionally  to  the  density  of  the  system,  it 
is  not  generally  used  in  flow  regimes  where  the  Navier-Stokes  equations  are  valid.  We 
will  demonstrate  quantitative  agreement  between  our  DSMC  simulations  and  analytic 
results  in  the  low-  and  high-pressure  limits. 

There  are  two  key  factors  which  make  DSMC  applicable  to  much  larger  systems 
than  those  tractable  by  MD.  First,  each  simulated  particle  can  represent  an  extremely 
large  number  of  gas  molecules.  In  the  calculations  performed  in  the  present  study,  the 
ratio  of  real  molecules  to  simulated  molecules  ranged  from  5.72  x  to  2.93  x  10^^. 
Second,  a  basic  assumption  in  DSMC  is  that  the  convective  motion  of  the  particles 
can  be  decoupled  from  the  intermolecular  collisions  on  a  timescale  comparable  to  the 
mean  collision  time.  This  results  in  a  much  longer  allowable  timestep  than  in  MD 
where  the  timestep  is  determined  by  the  intermolecular  force.  The  DSMC  code  used 
in  this  study  is  based  on  source  code  available  from  Bird,  with  modifications  made 
for  the  surface  reactions  described  below. 

B.  Gas-phase  model 

The  model  system  used  in  the  calculations  consists  of  a  gas  mixture  of  either 
H/H2  or  H/H2/CO  confined  in  one-dimension  between  a  diamond  growth  surface  at 
X  =  0  and  a  hot  activating  surface  at  x  =  L  as  shown  in  Fig.  1.  The  computational 
domain  is  divided  into  100  equally  spaced  cells.  The  growth  surface  is  maintained  at  a 
temperature  of  1200  K  for  all  cases  and  the  activating  surface  is  set  to  a  temperature 
in  the  range  2000-4000  K.  The  interactions  between  the  gas-phase  species  are  treated 
using  the  variable  soft  sphere  collision  model  [16,  17],  which  properly  accounts  for 
variations  in  the  viscosities  and  self-diffusion  constants  of  the  gas-phase  species  as  a 
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function  of  temperature.  Rotational  degrees  of  freedom  for  H2  and  CO  are  considered 
to  be  fully  equilibrated  at  the  local  gas  temperature.  A  discrete  level  model  for  the 
vibrational  state  of  H2  was  included,  but  its  effect  was  negligible  since  H2  was  found  to 
be  at  most  15%  vibrationally  excited  in  the  hottest  region  of  the  reactor.  Equilibration 
between  the  translational,  rotational,  and  vibrational  degrees  of  freedom  followed 
the  Larsen- Borgnakke  model  [18].  The  homogeneous  (gas-phase)  recombination  of 
hydrogen  atoms  was  originally  considered,  but  was  found  to  be  negligible  compared 
to  the  rate  of  recombination  on  the  diamond  surface.  The  parameters  for  the  gas- 
phase  species  are  listed  in  Table  1. 

C.  Surface  model 

The  accommodation  factor  for  all  surface  scattering  processes  is  set  to  be  unity, 
i.e.,  incident  molecules  become  completely  thermally  equilibrated  with  the  surface 
and  are  reflected  with  translational,  rotational,  and  vibrational  energies  chosen  from 
the  appropriate  Boltzmann  distribution.  While  this  assumption  is  usually  valid  for 
heavy  molecules,  light  species  such  as  H  and  H2  generally  have  accommodation  co¬ 
efficients  less  than  unity  because  of  the  large  mass  differences  between  the  incident 
molecules  and  the  surface  atoms.  A  large  majority  of  the  diamond  surface  sites  are 
hydrogen  terminated,  however.  Thus,  most  molecules  incident  on  the  surface  will  end 
up  striking  a  hydrogen  atom  rather  than  a  carbon  atom,  resulting  in  more  complete 
accommodation.  All  gas-surface  collisions  other  than  atomic  hydrogen  scattering  from 
the  diamond  growth  surface  and  molecular  hydrogen  scattering  from  the  activating 
surface  are  taken  to  be  non-reactive. 

The  growth  surface  consists  of  two  types  of  sites:  hydrogen-terminated  carbon 
sites  and  surface  radical  sites,  denoted  by  CdH  and  Cj,  respectively.  The  correspond¬ 
ing  fractional  coverages  of  the  growth  surface  by  the  two  types  of  sites  are  given 
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by  and  /*.  Two  elementary  gas-surface  reactions  involving  hydrogen  atoms  are 
considered:  the  abstraction  reaction 

CdH  +  H-.Q  +  H2,  (1) 

and  the  termination  reaction  [20] 

Q  +  H  ^  CdH.  (2) 

The  first  of  these  is  responsible  for  the  creation  of  surface  radical  sites  where  carbon 
addition  can  occur,  while  the  second  has  the  effect  of  passivating  the  growth  surface. 
Taken  together,  these  two  reactions  comprise  the  primary  mechanism  responsible  for 
the  conversion  of  atomic  to  molecular  hydrogen  under  CVD  conditions.  The  reverse 
processes  of  the  above  reactions  are  neglected.  Goodwin  [6]  has  estimated  on  the 
basis  of  thermochemistry  that  the  equilibrium  constant,  kffkr,  for  the  abstraction 
reaction  is  approximately  800,  while  Coltrin  and  Dandy  [21]  propose  a  surface  reaction 
probability  of  ~  4  x  10”“*  for  the  reverse  process.  The  thermal  desorption  of  hydrogen 
can  be  neglected  on  the  basis  of  the  large  activation-energy  barrier  (~75  kcal/mol 
[22])  for  the  process. 

Because  of  their  key  importance  in  diamond  growth,  the  rates  of  these  reactions 
have  been  of  considerable  interest.  Since  we  are  employing  a  particle-based  simulation 
technique,  it  is  convenient  to  express  the  gas-surface  rate  constants  in  terms  of  the 
surface  reaction  probabilities  jabs  and  jterm-  The  reaction  probability  can  be  related 
to  the  more  familiar  mass-action  kinetic  rate  constant  through 

where  T  is  the  surface-site  concentration  (3.14  x  10^®  sites  cm~^  for  diamond),  m 
is  the  sum  of  the  stoichiometric  coefficients  of  the  surface-phase  reactants,  and  W 
is  the  molecular  weight  of  the  gas-phase  species.  The  fractional  coverages  of  the 
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surface  can  be  expressed  in  terms  of  the  reaction  probabilities  through  the  expressions 
f  ~  ')terml (^Kterm  ~^'yabs')  s^ud  f  —  '^abs / {^term  ~^"fabs) •  The  probabilities  of  an  incident 
gas-phase  hydrogen  atom  abstracting  a  surface  hydrogen  or  terminating  a  surface 
radical  are  equal  and  given  by  'yterm'Jabs I {iterm  +  labs)- 

Experimentally  measured  rates  are  not  available  for  the  surface  reactions,  and  so 
estimates  from  analogous  gas-phase  reactions  or  numerical  simulations  must  be  used. 
The  results  presented  in  this  paper  were  obtained  using  the  values  of  'jabs  =  0.04 
and  'iterm  —  0.43  (at  1200  K)  proposed  by  Brenner  et  al.  [23],  based  on  molecular 
dynamics  simulations  employing  empirical  carbon- hydrogen  potentials  [24].  Other 
estimates  for  the  surface  reaction  probabilities  include  ')abs  =  0.22,  7term  =  1  (Harris 
[25]),  7oi,s  =  0.016,  7term  =  0.19  (Harris  and  Belton  [26]),  and  '^abs  =  0.0037,  7/erm 
=  0.016  (Frenklach  and  Wang  [27]).  The  parameters  from  the  molecular  dynamics 
simulations  were  chosen  because,  as  pointed  out  by  Goodwin  [6] ,  they  account  for  the 
important  physical  effects  of  both  energy  transfer  to  the  lattice  and  steric  hindrance. 

Although  the  rates  of  these  elementary  reactions  have  not  been  determined  ex¬ 
perimentally,  related  quantities  have  been  measured.  Harris  and  Weiner  [28]  cite  a 
value  of  0.12,  with  a  level  of  accuracy  of  about  a  factor  of  two,  for  the  probability 
of  hydrogen  atom  recombination  on  the  surface.  It  is  important  to  distinguish  be¬ 
tween  Jabs  and  7rec-  The  former  is  the  probability  that  a  hydrogen  atom  striking  a 
hydrogen-terminated  surface  site  participates  in  an  abstraction  reaction,  while  the 
later  is  the  likelihood  that  an  incident  hydrogen  atom  eventually  reacts  to  form  H2 
rather  than  scatter  non-reactively.  This  quantity  can  be  expressed  in  terms  of  the 
elementary  reaction  probabilities  as 

7»’ec  ~  ^'ytermKabs I (Siterm  “t"  7a6s)*  (4) 

Using  high-resolution  electron  energy  loss  spectroscopy  (HREELS),  Thoms  et  al.  [5] 
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were  able  to  determine  a  value  of  0.05  ±  0.01  for  the  ratio  'yabsllterm-  The  elementary 
reaction  probabilities  of  Brenner  et  al.  are  reasonably  consistent  with  both  of  these 
experimental  studies. 

The  diamond  surface  chemistry  described  above  results  in  the  conversion  of  atomic 
to  molecular  hydrogen.  For  diamond  growth  to  occur,  a  constant  flux  of  hydrogen 
atoms  to  the  growth  surface  is  required.  To  provide  a  source  of  atomic  hydrogen,  H2 
molecules  incident  on  the  activating  surface  are  allowed  to  dissociate  with  a  probabil¬ 
ity  Idiss-  Most  of  the  results  presented  in  this  paper  were  obtained  with  set  equal 
to  unity.  The  degree  of  hydrogen  dissociation  under  CVD  conditions  is  not  always 
known  and  this  choice  for  '^diss  allows  us  to  establish  an  upper  limit  on  the  H  atom 
flux  for  a  given  set  of  operating  conditions. 

III.  NUMERICAL  RESULTS 

In  this  section,  we  present  the  results  of  our  DSMC  calculations.  Since  the  pa¬ 
rameter  space  of  interest  is  too  large  to  study  completely,  we  have  chosen  a  set  of 
parameter  values  as  standard  conditions  and  looked  at  the  effect  of  varying  one  pa¬ 
rameter  at  a  time  from  this  set.  Our  standard  conditions  comprise  a  gas  mixture  of 
molecular  and  atomic  hydrogen  with  L  =  1  cm,  Tsuh  =  1200  K,  Tact  =  3000  K,  'fterm 
=  0.43,  Jabs  =  0.04,  and  jdiss  =  1-0.  The  effect  of  varying  each  parameter  was  inves¬ 
tigated  over  a  range  of  reactor  pressures  P.  This  allows  us  to  compare  the  numerical 
results  to  the  analytic  expressions  in  both  the  low-  and  high-pressure  limits.  We  can 
also  determine  the  range  of  validity  of  the  analytic  forms  and  obtain  quantitative 
results  in  the  transition  regime. 

In  all  of  the  calculations,  the  system  is  initialized  with  a  linear  temperature  gra¬ 
dient  extending  across  the  length  of  the  computational  domain.  The  gas-phase  com¬ 
position  is  initially  uniform  everywhere  with  atomic  and  molecular  hydrogen  mole 
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fractions  of  2/3  and  1/3,  respectively.  The  temperature  and  species  gradients  are  al¬ 
lowed  to  evolve  until  the  steady-state  conditions  are  reached.  Statistics  are  collected 
over  20,000-80,000  DSMC  timesteps  only  after  it  is  clear  that  the  temperature  and 
concentration  gradients  are  no  longer  evolving.  The  pressure  of  the  system  is  not 
set  since  we  cannot  predict  a  priori  the  total  number  of  particles  in  the  final  state. 
Instead,  we  set  the  initial  gas  density  and  determine  the  pressure  at  the  end  of  the 
calculation  by  evaluating  the  net  momentum  transfer  to  the  bounding  surfaces.  In 
all  cases,  the  pressures  measured  at  the  two  walls  agreed  to  within  less  than  0.1%. 
The  difference  between  the  heat  fluxes  at  the  substrate  and  activating  surface  ranged 
from  0.1%  at  low  pressures  to  at  most  1%  at  high  pressures. 

The  main  focus  of  this  work  is  the  effect  of  the  reactor  operating  conditions  on 
the  hydrogen  atom  flux  (^h)  to  the  growth  surface.  Also  of  interest  is  the  partition¬ 
ing  of  the  total  heat  flux  (qtot)  to  the  substrate  into  contributions  due  to  thermal 
conduction  (qcond)  and  hydrogen  atom  recombination  (qrec)-  The  recombinative  heat 
flux  is  related  to  the  H-atom  flux  through  the  expression 

qrec  ~  4^V{'^Tec^H^ (5) 

where  =  104.2  kcal  mol“^  is  the  heat  of  formation  of  atomic  hydrogen  from 

molecular  hydrogen.  The  conductive  heat  flux  qcond  is  defined  as  the  difference  be¬ 
tween  the  total  (kinetic  -f-  rotational  -f  vibrational)  energies  of  the  incident  and  re¬ 
flected  molecules.  The  energy  flux  due  to  the  conversion  of  two  hydrogen  atoms  with 
six  degrees  of  freedom  into  a  single  hydrogen  molecule  with  five  degrees  of  freedom  is 
accounted  for  by  qcond- 


A.  Effect  of  separation  between  substrate  and  heat  source 
In  Fig.  2,  the  species  and  temperature  profiles  for  an  H/H2  mixture  at  pressures  of 
0.050,  0.45,  and  6.6  Torr  are  shown.  At  low  pressures,  rarefied  gas  effects  are  clearly 
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observed.  Most  notably,  the  temperature  profile  does  not  vary  continuously  from  1200 
K  to  3000  K  as  predicted  from  continuum  theory,  but  exhibits  jump  discontinuities  at 
the  surfaces  typical  of  high  Kn  flows.  The  temperature  jump  is  a  direct  consequence 
of  the  extremely  low  density  of  the  gas,  since  individual  gas  molecules  travel  relatively 
large  distances  before  undergoing  thermalizing  collisions.  As  the  pressure  is  increased, 
the  continuum  regime  is  approached  and  the  temperature  jump  is  diminished.  The 
effect  of  pressure  on  the  species  profiles  is  also  quite  profound.  At  0.050  Torr,  the 
hydrogen  is  almost  completely  dissociated  and  there  is  little  spatial  variation  in  the 
gas  phase  composition.  As  the  pressure  is  increased,  the  counter-diffusion  of  the  two 
species  is  retarded  and  larger  concentration  gradients  develop. 

The  hydrogen-atom  flux,  total  heat  flux,  and  the  contributions  from  thermal 
conduction  and  hydrogen  recombination  are  shown  in  Fig.  3  as  a  function  of  pressure 
for  values  of  L  equal  to  0.25,  0.5,  and  1.0  cm.  In  the  low-pressure  limit,  qcond  and 
qrec  exhibit  a  linear  dependence  on  pressure  and  are  independent  of  L.  At  higher 
pressures,  the  components  of  the  heat  flux  asymptote  to  constant  values  that  are 
strongly  dependent  on  L.  At  a  pressure  in  the  range  of  0. 5-2.0  Torr,  a  crossover 
between  qcond  and  qrec  occurs.  In  the  low-pressure  limit,  qcond  is  greater  than  qreci 
while  at  25  Torr,  qrec  exceeds  qcond  by  a  factor  of  2-4. 

B.  Effect  of  CO  dilution 

In  Fig.  4,  the  species  and  temperature  profiles  for  the  H/H2/CO  mixture  with  a 
CO-H  stoichiometry  of  1:4  are  shown  for  three  different  pressures.  The  temperature 
and  hydrogen  profiles  show  the  same  trends  as  were  seen  in  the  H/H2  mixture,  but 
small  quantitative  differences  exist  between  the  temperature  profiles.  A  comparison 
of  Figs.  2a  and  4a  shows  that  the  temperature  jump  at  the  walls  is  slightly  less  for 
the  H/H2/CO  mixture.  This  is  due  to  the  fact  that  the  addition  of  carbon  monoxide 
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to  the  gas  phase  reduces  the  average  molecular  mean-free  path.  The  influence  of 
thermal  diffusion  (Soret  effect),  the  preferential  diffusion  of  lighter  species  toward 
regions  of  higher  temperature,  can  be  seen  in  the  CO  profiles.  The  CO  mole  fraction 
for  the  simulation  carried  out  at  6.5  Torr,  e.g.,  varies  from  0.3  at  the  substrate  to 
0.2  at  the  activating  surface.  Additional  calculations  were  performed  in  which  the 
mass  of  CO  was  reduced  to  that  of  atomic  hydrogen,  or  the  surface  reactions  were 
turned  off  to  eliminate  the  counter- flows  of  H  and  H2.  The  results  confirm  that  the 
CO  concentration  gradient  is  due  to  thermal  diffusion  driven  by  the  disparate  masses 
of  CO  and  H/H2. 

The  hydrogen  atom  and  thermal  fluxes  are  shown  in  Fig.  5  for  the  H/H2/CO 
mixture.  Also  displayed  for  comparison  are  the  results  for  the  H/H2  mixture.  The 
general  trends  and  behavior  of  qtou  9rec?  and  qcond  are  the  same  for  the  two  cases, 
with  both  qcond  and  g^ec  consistently  lower  for  the  H/H2/CO  mixture.  Although  CO 
comprises  only  about  25%  of  the  gas  mixture  by  mole  fraction,  the  effects  on  the  fluxes 
are  much  larger.  The  high-pressure  limit  of  qcond  is  reduced  by  roughly  a  factor  of  two. 
The  recombinative  heat  flux  has  not  yet  reached  its  asymptotic  value,  but  it  appears 
to  be  reduced  by  a  factor  in  the  range  1.5-2.  Because  of  its  large  mass  and  size, 
as  compared  to  atomic  and  molecular  hydrogen,  the  carbon  monoxide  significantly 
retards  the  counter-diffusion  of  the  hydrogen  species.  This  leads  to  a  severe  reduction 
of  both  qcond  ^nd  qrec  * 

C.  Effect  of  temperature  at  activating  surface 

A  series  of  calculations  were  carried  out  to  determine  the  effect  of  variations  in 
the  temperature  at  the  activating  surface  on  the  hydrogen-atom  and  heat  fluxes.  The 
results  are  shown  in  Fig.  6  for  the  H/H2  mixture  with  L  =  I  cm.  As  expected 
from  continuum  theory  assuming  a  dependence  of  the  thermal  conductivity,  the 
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^  c /3 

conductive  heat  flux  in  the  high-pressure  limit  is  roughly  proportional  to  {T^ct  ” 
^sub)'  high  pressures,  qcond  rises  steadily  but  more  slowly  with  Tact-  9rec  likewise 
increases  monotonically  with  Tact  high  pressures,  but  the  changes  are  relatively 
small.  Over  the  temperature  range  2000  ^  Tact  <  4000,  qrec  varies  by  only  about  25%. 
In  contrast,  at  low  pressures  qrec  shows  the  counterintuitive  behavior  of  decreasing  as 
Tact  is  increased.  This  is  due  to  the  fact  that  the  increased  thermal  velocities  of  the 
molecules  are  more  than  offset  by  the  lower  gas  density  as  will  be  shown  in  the  next 
section.  For  our  choice  of  operating  conditions,  the  dependence  of  qrec  on  Tact  crosses 
over  from  the  low-  to  high-pressure  behaviors  for  pressures  in  the  range  of  a  few  Torr. 

D.  Effect  of  H2  dissociation  probability 

All  of  the  results  presented  to  this  point  were  obtained  assuming  a  100%  dis¬ 
sociation  of  molecular  hydrogen  at  the  activating  surface.  The  various  activation 
techniques  used  in  CVD  reactors  can  lead  to  the  direct  dissociation  of  H2  through 
mechanisms  other  than  neutral  gas-phase  chemistry  -  e.g.,  electron  impact  dissocia¬ 
tion  in  a  plasma  -  resulting  in  elevated  H/H2  ratios.  Nonetheless,  100%  dissociation 
of  H2  is  clearly  an  optimal  upper  bound.  In  this  section  we  consider  the  effect  of 
incomplete  dissociation  of  H2  at  the  activating  surface.  Figure  7  compares  the  results 
for  H2-H  mixtures  with  L  equal  to  1  cm  for  'ydiss  set  to  1.0,  0.5,  and  0.1. 

The  conductive  component  of  the  heat  flux  is  very  insensitive  to  the  value  of 
'ydiss-  Reducing  ydiss  from  1.0  to  0.5  results  in  negligible  changes  in  qcond^  while 
further  decreasing  ydiss  to  0.1  leads  to  at  most  a  few  percent  reduction  in  qcond-  A 
detailed  explanation  for  this  behavior  will  be  given  in  §IV,  but  it  can  be  at  least 
be  partially  understood  in  terms  of  the  species  profiles  for  the  ydiss  =  1.0  and  0.1 
cases  shown  in  Fig.  8.  In  the  low-pressure  limit,  the  degree  of  hydrogen  dissociation 
depends  strongly  on  ydiss-  Because  the  molar  heat  fluxes  for  atomic  and  molecular 
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hydrogen  in  the  free-streaming  limit  differ  by  less  than  20%,  the  effect  on  the  thermal 
conductivity  of  the  mixture  is  relatively  small.  At  high  pressures,  the  composition 
and  hence  the  thermal  conductivity  of  the  gas  depends  only  weakly  on  'ydiss-  Although 
qrec  shows  a  stronger  dependence  on  'jdiss  than  does  qcondi  the  variation  of  qrec  with 
'jdiss  is  still  rather  small  until  •jdiss  is  reduced  to  0.1. 

E.  Effect  of  surface  reaction  probabilities 

It  was  pointed  out  in  §IIc  that  there  is  considerable  uncertainty  in  the  probabil¬ 
ities  of  the  gas-surface  reactions.  None  of  the  reaction  rates  are  at  present  available 
experimentally,  so  the  rates  must  be  either  derived  from  first  principles  or  taken  from 
analogous  gas  phase  reactions.  While  the  values  for  'jabs  and  'yterm  proposed  by  various 
investigators  vary  by  almost  two  orders  of  magnitude,  the  ratio  ')termllabs  lies  within 
a  much  narrower  range.  The  values  proposed  by  Brenner  et  al.  are  consistent  with 
the  hydrogen  recombination  probability  7rec  measured  by  Harris  and  Weiner  [28],  but 
this  agreement  does  not  guarantee  the  validity  of  their  values. 

In  Fig.  9,  we  compare  the  results  of  simulations  on  the  H/H2  mixture  for  three 
different  sets  of  surface  reaction  probabilities.  In  the  first  case,  our  standard  values 
from  Brenner  et  al.  [23,  24]  {■yterm  =  0.43  and  'yaks  =  0.04)  are  used.  In  the  second 
case,  the  probabilities  of  both  surface  reactions  are  doubled.  Since  the  fractions  of  the 
surface  sites  in  the  two  types  of  configurations,  CdH  and  Cj,  are  determined  solely 
by  the  ratio  'ytermHabs,  changing  the  surface  reaction  probabilities  in  this  way  has  no 
effect  on  the  state  of  the  surface.  Finally,  in  the  third  case,  -fabs  is  doubled  while  'yterm 
is  left  unchanged. 

Variations  in  the  surface  reaction  probabilities  have  virtually  no  effect  on  qcond- 

♦ 

The  explanation  for  this  behavior  is  identical  to  that  given  in  the  previous  section.  In 
the  low-pressure  limit,  the  heat  flux  depends  very  weakly  on  the  degree  of  hydrogen 
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dissociation,  while  in  the  high-pressure  limit,  the  gas-phase  composition  is  nearly 
independent  of  the  surface  reaction  probabilities.  The  quantities  qrec  and  (l>i\  depend 
strongly  on  jabs,  but  are  nearly  independent  of  the  value  for  jterm-  This  is  because 
the  total  recombination  probability,  jrec  =  ‘^Jabsltermfilabs  +  Jterm),  is  approximately 
^Jabs  so  long  as  7term  ^  labs-  At  high  pressures,  the  recombinative  heat  fluxes  for  the 
three  cases  appear  to  be  converging,  while  at  low  pressures  the  hydrogen  atom  flux 
is  independent  of  jabs- 


IV.  ANALYSIS 

Our  model  system  can  be  analyzed  to  derive  expressions  for  the  particle  and  en¬ 
ergy  fluxes  in  both  the  low-  and  high-pressure  regimes.  In  this  section,  these  results 
are  presented  and  then  used  to  interpret  and  understand  further  the  DSMC  simula¬ 
tions  described  above.  Details  of  the  derivations  are  relegated  to  the  Appendix. 


A.  Low-pressure  regime 

At  low  pressures,  the  gas  can  be  treated  as  a  mixture  of  counter-streaming  pop¬ 
ulations  which  collide  only  with  the  bounding  surfaces.  Requiring  the  oppositely 
directed  mass  fluxes  to  balance  at  those  surfaces  fixes  the  relative  populations  of  cold 
and  hot  particles  of  each  species  separately  and  those  of  H  atoms  and  H2  molecules 
jointly.  The  results  are 


4>u  =  n{2kBlTrm}if^‘^ 


rpl/2nnl/2 

act  sub 


P  =  nkeT^iX!, 


Xa, 


Ml  “  7ree) 


1/2 


Xh 


+ 


Tll^  +  {I  -  ldiss)TlS 

qcond  =  {2li:fl^nkTTliXil  '  TI!2) 


JVhj  +  A’co 


(6) 


(7) 


14 


X  {2m^''^XY{  +  +  Zmco^Xco) 

+  <i>HlreckB  {Tact  “  Tal^T^l^  +  Tsub^  /2, 


Qrec  —  4^W'^rec^BTbindl‘^i 


Tree  _  ol/2  Idiss  V-  /inN 

fi-'Y  ')ri/2  .  ji/2  T'/'  +  ri -'v^-  'IT'/' 

In  these  equations,  and  Xs  are  the  mass  and  fractional  concentration,  respectively, 
of  species  s,  and  T^ind  =  XH^/ks  =  5.25  x  10“*  K  is  the  binding  energy  of  the 
H2  molecule  expressed  as  a  temperature.  The  contributions  of  H2  and  CO  to  the 
conductive  heat  flux  in  (8)  are  enhanced  over  that  of  H  by  a  factor  of  3/2,  due 
to  excitation  of  the  two  rotational  degrees  of  freedom  of  the  diatomic  molecules. 
Excitation  of  the  vibrational  degrees  of  freedom  of  H2  and  CO  was  ignored  in  the 
analysis,  in  light  of  our  DSMC  results.  The  recombination  and  dissociation  reactions 
contribute  through  the  conversion  at  the  substrate  of  two  H  atoms  with  six  degrees  of 
freedom  into  a  single  H2  molecule  with  five  degrees  of  freedom,  and  the  inverse  process 
at  the  activating  surface.  If  these  reactions  do  not  occur,  i.e.,  7rec  =  'ydiss  =  0,  then 
the  expressions  above  reduce  to  those  given  elsewhere  [3],  or  to  trivial  identities. 

The  particle  and  energy  fluxes  and  the  pressure  are  all  linear  in  the  number 
density  in  the  free-streaming  regime.  Consequently,  if  the  temperatures  and  reaction 
probabilities  are  held  fixed,  the  fluxes  will  scale  linearly  with  the  pressure.  The  fluxes 
also  are  predicted  to  be  independent  of  the  separation  between  the  substrate  and  the 
activation  source.  Both  of  these  trends  are  consistent  with  the  low-pressure  behavior 
of  our  DSMC  simulations.  A  comparison  between  the  lowest  pressure  H/H2  DSMC 
simulations  and  the  analytic  forms  is  given  in  Table  2.  The  quantitative  agreement 
between  the  two  approaches  is  very  good.  It  is  worth  noting  that  the  discrepancy 
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in  the  value  for  qcond  is  smallest  for  the  largest  Knudsen  number  case,  L  =  0.25  cm. 
As  the  separation  between  the  substrate  and  the  activating  surface  is  increased,  with 
pressure  held  fixed,  collisional  effects  become  more  important  and  the  free-streaming 
approximation  loses  validity.  The  effect  of  L  on  the  transition  from  a  collisionless  to 
a  collisional  regime  can  also  be  seen  in  Fig.  3b.  As  L  is  decreased,  the  fluxes  exhibit 
linear  scaling  over  a  larger  pressure  range. 

The  large  dissociation  fraction  in  our  standard  case  is  due  to  the  high  probability 
of  H2  dissociation  at  the  activating  wall  relative  to  that  of  H  atom  recombination  at 
the  substrate.  Rearranging  Eq.  (10)  yields 


V  _  0-1/2  7rec  +  (1  -  'ydiss)  Till  V 


(11) 


Thus,  ')rec  •<  Idiss  <  1  implies  C  Xu,  and  quite  generally  the  relative  concen¬ 
trations  are  very  sensitive  to  the  relative  surface  reaction  probabilities.  We  observed 
these  trends  in  our  calculations  with  variable  reaction  rates,  particularly  those  shown 
in  Fig.  8.  At  the  same  time,  we  found  the  conductive  heat  fluxes  in  Figs.  7  and  9  to  be 
essentially  independent  of  the  H/H2  ratio.  This  is  because  the  thermal  conductivities 
of  H  and  H2,  as  inferred  from  Eq.  (8),  are  nearly  identical:  the  larger  mass  of  the 
molecule  is  compensated  by  its  extra,  rotational  degrees  of  freedom.  The  recombi- 
native  heat  flux,  in  contrast,  shows  two  distinct  trends.  As  shown  in  Fig.  7,  it  is 
sensitive  to  the  dissociation  probability  'ydUs  only  when  the  latter  is  very  small.  This 
follows  immediately  from  Eq.  (11),  since  for  jdiss  ~  7rec,  the  H2  and  H  concentrations 
become  comparable.  The  resulting  drop  in  Xu  produces  a  proportional  decline  in  the 
hydrogen-atom  and  recombinative  heat  fluxes.  So  long  as  7rec  “ydiss,  the  hydrogen 
concentration  is  insensitive  to  these  probabilities.  The  recombinative  heat  flux,  how¬ 
ever,  is  directly  proportional  to  q^ec,  as  indicated  by  Eq.  (9)  and  illustrated  in  Fig.  9. 
Both  of  these  trends  just  reflect  the  fact  that  the  slow  recombination  of  H  atoms  at 
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the  substrate  is  a  rate-limiting  step  in  this  process. 

We  noted  in  Fig.  6  that  varying  the  temperature  of  the  activating  surface  has 
a  pronounced  effect  on  the  conductive  heat  flux  but  a  very  limited  impact  on  the 
hydrogen  atom  and  recombinative  heat  fluxes.  As  a  first  approximation,  we  can  reduce 
Eqs.  (7)  and  (8)  to  />  =  nksT^'?T',S  and  -  tUI), 

respectively.  This  leads  to  the  result 

qcond^P{Tli^-Tlil).  (12) 

Using  this  expression  we  find  that  the  ratio  qcond{Tact  =  4000) /^coni (Tact  =  2000)  is 
equal  to  2.8,  in  excellent  agreement  with  the  lowest  pressure  numerical  results  which 
yield  a  ratio  of  2.7.  Retaining  the  terms  that  have  been  omitted  in  the  derivation 
of  Eq.  (12)  further  improves  the  agreement  between  the  two  sets  of  results.  These 
smaller  terms  are  important  in  the  calculation  of  the  recombinative  heat  flux,  whose 
temperature  dependence  is  more  subtle  and  requires  a  more  careful  derivation.  Using 
Eqs.  (6),  (7),  (9),  and  (10)  we  find 

(13) 


Although  the  effect  is  relatively  small,  increasing  Tact  for  fixed  pressure  leads  to  a 
decrease  in  qrec  as  shown  in  Fig.  6b.  The  increase  in  the  average  thermal  speed  with 
increasing  Tact  in  Eq.  (6)  is  more  than  offset  by  the  simultaneous  reduction  in  the 
total  number  density. 

Under  our  standard  conditions  it  happens  that  the  conductive  heat  flux  is  greater 
than  the  recombinative  component.  Making  the  reasonable  approximations  q^ec  •C  1, 
Xhj  <C  A’h,  and  T^ind  ^  Tact,  we  find  from  Eqs.  (8)  and  (9)  that 


^cond 

Qrec 


Tact  —  Taub 
'yrecThindl^ 


(14) 
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Thus,  as  is  illustrated  in  Figs.  6  and  9,  the  flux  ratio  is  sensitive  to  the  temperature 
of  the  activating  surface  (through  its  effect  on  the  magnitude  of  the  conductive  flux) 
and  to  the  recombination  probability  of  H  atoms  at  the  substrate  (which  affects  the 
size  of  the  recombinative  flux). 

Finally,  adding  carbon  monoxide  as  a  diluent  at  low  pressures  has  the  principal 
effect  of  reducing  both  the  atomic  and  molecular  hydrogen  concentrations  propor¬ 
tionately,  and  therefore  the  H  atom  flux  and  the  heat  fluxes.  CO  contributes  to  the 
conductive  flux  much  more  weakly  due  to  its  larger  mass  (and  consequently  longer 
transit  time  between  the  activating  surface  and  substrate).  Thus,  a  reduction  in  the 
heat  fluxes  of  between  12.5%  and  25%  -  depending  upon  the  extent  to  which  H2  is 
dissociated  in  the  gas  -  is  expected.  The  observed  effects  in  Fig.  5  are  somewhat 
larger  than  this,  evidently  due  to  the  effects  of  rare  but  not  inconsequential  collisions 
with  CO  molecules.  We  noted  the  latter  in  conjunction  with  the  temperature  profiles 
of  Fig.  4. 

B.  High-pressure  regime 

At  high  pressures,  the  gas  can  be  analyzed  using  a  hydrodynamical  description. 
This  generally  involves  the  solution  of  a  set  of  coupled,  highly  nonlinear  partial  differ¬ 
ential  equations.  As  detailed  in  the  Appendix,  a  few  simplifying  assumptions  permit 
us  to  draw  some  analytically  tractable  conclusions  which  turn  out  to  be  quantitatively 
useful.  The  approximate  results  for  a  pure  H/H2  mixture  are 

^  (15) 

qcond  «  {Tact  -  Tsub)  I^L  +  Z4>HlreckB  {Tact  +  T,ub)  /8,  (16) 

Qrec  ~  ^HKrecksTbind (1"^) 
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and 


rp-\l2  ^^5^6  rp-\l2  y^act 

IrecJ-suh  IdissJ- act 

?»  (2xmHfcs)^'^^ln2  (c?h-H2  +  ^h-Hj) 

In  these  equations,  «?h-H2  =  ^^h-H2  ?  where  Dh-H2  is  the  coefficient  for  binary  diffu¬ 
sion  of  H  atoms  through  H2  molecules,  and  kh  ^•nd  /chj  ars  the  thermal  conductivi¬ 
ties  of  atomic  and  molecular  hydrogen.  All  of  these  transport  coefficients  vary  with 
temperature,  which  has  been  taken  into  account  in  an  approximate  way  by  simply 
averaging  the  values  at  the  substrate  and  activating  surface.  For  our  standard  case, 
the  average  diffusion  constant  and  thermal  conductivity  are  2.0  x  10^°  cm“^  s“^  and 
1.2  X  10“^  W  cm“^  K“^,  respectively. 

The  fluxes  (15)-(17)  are  independent  of  pressure,  indicating  that  a  saturation  of 
the  fluxes  occurs  as  the  pressure  increases  from  the  free-streaming  regime.  They  also 
depend  inversely  on  the  separation  L  between  the  substrate  and  activating  surface, 
the  only  macroscopic  scale  length  in  our  system.  These  are  very  general  features  of  the 
hydrodynamic  description,  independent  of  any  other  simplifying  approximations  that 
have  been  made  to  derive  the  above  expressions.  The  two  behaviors  are  exhibited  by 
our  DSMC  results  in  the  high-pressure  limit,  particularly  those  for  the  conductive  heat 
fluxes.  The  hydrogen-atom  and  recombinative  heat  fluxes  have  not  quite  saturated  at 
the  maximum  pressures  assumed  in  the  DSMC  runs,  but  the  trend  in  this  direction 
is  fairly  clear.  Quantitatively,  we  compare  results  from  our  highest-pressure  DSMC 
simulations  with  the  analytic  forms  in  Table  3.  Given  the  approximations  inherent 
in  our  analysis,  the  agreement  with  the  DSMC  values  is  very  good.  As  is  to  be 
expected,  the  closest  correspondence  is  obtained  for  the  largest  gap,  and  thus  the 
most  collisional  system,  L  =  I  cm. 

All  of  the  fluxes  depend  implicitly  on  the  temperature  of  the  activating  surface 
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through  the  transport  coefficients  for  thermal  conduction  and  binary  diffusion.  These 
coefficients  rise  monotonically  over  the  temperature  range  of  interest.  Thus,  as  is 
evident  in  Fig.  6,  the  atomic  hydrogen  and  recombinative  heat  fluxes  increase  as  the 
activating  surface  temperature  increases.  The  numerical  values  rise  by  about  25%  as 
Tact  varies  from  2000  to  4000  K.  From  the  atomic  data,  an  increase  of  about  40%  in  the 
diffusion  coefficient,  and  thus  the  fluxes,  is  expected.  This  discrepancy  probably  just 
reflects  the  clearly  incomplete  convergence  of  the  DSMC  fluxes  to  their  asymptotic 
values.  The  conductive  heat  fluxes,  which  have  more  nearly  converged,  increase  due 
to  both  the  greater  temperature  difference  between  substrate  and  activating  surface 
and  the  enhanced  conductivity.  Here,  we  can  analytically  estimate  8.2  and  40.0  W 
cm“^  for  qcondiTact  =  2000 A")  and  qcond{Tact  =  4000 A"),  respectively,  in  quite  good 
agreement  with  the  values  in  Fig.  6. 

The  analysis  predicts  that  the  heat  and  H-atom  fluxes  are  independent  of  the 
dissociation  probability  for  H2  molecules  at  the  activating  surface.  As  seen  in  Fig.  7, 
this  is  true  to  a  high  degree  of  accuracy  for  the  conductive  heat  fluxes  and  is  true  of 
the  H-atom  and  recombinative  heat  fluxes  except  at  very  low  dissociation  probability, 
7<itss  =  0.1.  Careful  inspection  of  the  figure  suggests,  however,  that  all  of  the  curves 
for  these  fluxes  are  indeed  converging  as  the  pressure  increases.  Indeed,  it  appears 
that  the  divergence  of  the  curves  just  reflects  a  slower  approach  to  the  asymptotic 
fluxes  when  'ydiss  is  small.  Specifically,  we  note  from  (18)  that  the  H2  concentration 
at  the  activating  surface  declines  to  its  asymptotic  value  of  zero  as  {'ydissP)~^-  Thus, 
the  convergence  with  increasing  pressure  of  the  species  profiles,  and  consequently  the 
particle  fluxes,  is  slower  for  smaller  values  of  the  dissociation  probability.  The  same 
statement  is  true  with  respect  to  the  recombination  probability  of  H  atoms  at  the 
substrate.  Evidently,  it  is  in  fact  the  very  small  value  of  7rec  (=  0.073)  that  delays 
the  convergence  of  the  species  profiles  and  the  saturation  of  the  fluxes  in  all  of  the 
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cases  we  have  examined.  As  corroboration  for  this  assertion  we  refer  to  the  species 
profiles  in  Fig.  8,  which  show  that  the  H  concentration  at  the  substrate  is  still  about 
20%  at  24  Torr  of  pressure,  while  the  H2  concentration  at  the  activating  surface  is 
nearly  zero. 

Varying  the  recombination  probability  7rec  has  an  effect  on  the  heat  fluxes  simi¬ 
lar  to  that  of  varying  'ydiss*  doubling  ^rec  has  no  discernible  effect  on  the  conductive 
heat  fluxes,  and  seems  only  to  hasten  the  saturation  of  the  recombinative  flux  as  the 
pressure  increases,  without  changing  its  asymptotic  value.  The  H-atom  flux  to  the 
substrate,  on  the  other  hand,  is  reduced  by  a  factor  of  two  (see  Fig.  9).  This  result  is 
in  accord  with  the  analytical  result  (15),  wherein  is  inversely  proportional  to  7rec* 
Physically,  this  happens  because  the  rate  at  which  H  atoms  can  be  transported  to  the 
substrate  is  limited  by  the  binary  diffusion  coefficient;  increasing  the  recombination 
probability  simply  reduces  the  (asymptotically  small)  concentration  of  H  atoms  at 
the  substrate  until  the  recombination  and  diffusion  rates  come  into  balance.  It  is 


interesting  to  note  that  the  effect  of  7rec  at  high  pressures  is  the  opposite  of  that  at 
low  pressures,  where  it  had  no  effect  on  the  H-atom  flux  but  contributed  proportion¬ 
ately  to  the  associated  recombinative  flux.  This  is  reflected  also  in  the  ratio  of  the 


conductive  and  recombinative  heat  fluxes. 


qcond  ^  ~  '^sub 


qrec  ^H-H2  ^  ksTund 

whose  independence  of  7rec  contrasts  sharply  with  the  flux  ratio  (14)  at  low  pressures. 


V.  SUMMARY  AND  CONCLUSIONS 

One-dimensional  DSMC  calculations  have  been  carried  out  for  the  purpose  of 
studying  diffusion-dominated  diamond  CVD  processes.  The  effects  of  pressure,  dis¬ 
tance  between  the  substrate  and  the  activating  surface,  gas-phase  composition,  re¬ 
actor  temperature,  surface  reaction  probabilities,  and  degree  of  molecular  hydrogen 
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dissociation  on  the  hydrogen- atom  and  thermal  fluxes  to  the  diamond  surface  were 
investigated.  To  supplement  and  aid  in  the  interpretation  of  the  numerical  results, 
analytic  expressions  were  derived  for  the  low-  and  high-pressure  limits.  Good  quan¬ 
titative  and  qualitative  agreement  was  found  between  the  two  approaches. 

Several  of  the  model  parameters  were  shown  to  play  a  critical  role  in  determining 
the  heat  and  hydrogen-atom  fluxes  to  the  surface.  In  many  instances,  the  effect  of 
varying  a  particular  parameter  depends  sensitively  upon  the  pressure  regime  in  which 
one  operates.  At  low  pressures,  the  H-atom  and  heat  fluxes  all  scale  directly  with 
the  gas  pressure,  ?.e.,  the  number  density  of  particles.  The  greater  the  number  of 
carriers  (of  mass,  momentum,  or  energy),  the  larger  the  fluxes  that  can  be  supported 
by  streaming  of  particles  between  the  activating  surface  and  the  substrate.  At  high 
pressures,  in  contrast,  the  fluxes  eventually  saturate  at  fixed  values.  In  this  collisional 
regime,  increasing  the  number  of  carriers  by  raising  the  pressure  is  exactly  compen¬ 
sated  for  by  the  concomitant  increase  in  the  density  of  scatterers,  which  impedes  the 
transport  of  mass,  momentum,  and  energy  across  the  gap. 

A  reciprocal  role  is  played  by  the  distance  separating  the  activating  surface  and 
the  substrate.  The  particles  free-stream  between  the  two  at  very  low  pressures,  so 
the  fluxes  are  independent  of  the  separation.  At  high  pressures,  on  the  other  hand, 
the  fluxes  are  driven  by  the  gradients  in  the  temperature  and  composition  across  the 
gap.  The  strength  of  these  gradients  and  their  associated  thermodynamic  forces  vary 
inversely  as  the  width  of  the  gap. 

We  examined  the  effect  of  oxygen-bearing  species  in  diamond  CVD  by  adding  car¬ 
bon  monoxide  to  the  hydrogen  mixture.  CO  reduces  all  of  the  fluxes  to  the  substrate 
and  activating  surface.  In  the  low-pressure  limit,  the  total  fluxes  are  just  a  linear 
superposition  of  independent  contributions  from  each  species.  As  a  result,  the  H  and 
H2  fluxes  are  diminished  exactly  in  proportion  to  the  CO  mole  fraction.  At  higher 
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pressures,  CO  has  the  further  effect  of  substantially  retarding  the  counter-diffusion 
of  atomic  and  molecular  hydrogen  through  collisions  with  those  species.  Thus,  the 
detrimental  influence  of  CO  addition  on  the  H-atom  transport  rises  as  the  pressure 
increases. 

Varying  the  temperature  of  the  activating  surface  produced  only  subtle  effects 
on  the  hydrogen-atom  and  recombinative  heat  fluxes.  At  low  pressures,  these  fluxes 
exhibit  the  counter-intuitive  behavior  of  decreasing  as  the  temperature  Tact  increases 
at  fixed  pressure.  This  is  due  to  the  decline  in  the  number  density  of  the  gas  as 
the  average  temperature  is  raised.  A  cross-over  occurs  at  intermediate  pressures,  to 
the  high-pressure  behavior  of  the  fluxes  rising  slowly  with  rising  Tact-  The  increased 
efficacy  of  diffusion  of  H  through  H2  at  higher  temperatures  is  responsible  for  this 
(weak)  trend. 

As  one  would  expect,  in  contrast  the  activating  temperature  has  a  very  strong 
influence  on  the  conductive  heat  flux,  qcond  increases  monotonically  with  increasing 
Tact,  from  the  rather  slow  dependence  as  {Talt  ~  TIH)  in  the  free-streaming  regime 
at  low  pressures,  to  the  much  steeper  dependence  as  roughly  {T^if  —  Tilt)  in  the  con¬ 
tinuum  regime  at  high  pressures.  The  latter  reflects  an  approximate  T'^!^  dependence 
of  the  thermal  conductivity  of  the  gas  on  its  temperature. 

In  addition  to  the  operating  parameters  which  can  be  varied  and  monitored  in  an 
experimental  CVD  reactor,  we  also  have  studied  the  role  of  the  surface  reaction  and 
dissociation  rates.  The  probabilities  of  H-atoms  terminating  a  reactive  site,  '^term, 
and  abstracting  another  H  atom  from  a  previously  terminated  site  to  form  H2,  'labs, 
jointly  determine  both  the  state  of  the  surface  and  the  total  probability  of  H-atom 
recombination  on  the  surface,  7rec  =  2'yabs'ftermli'iabs  +  7term)'  Although  the  absolute 
values  of  these  reaction  probabilities  are  quite  uncertain,  there  is  a  consensus  that 
the  termination  reaction  is  much  faster  than  the  abstraction  reaction,  '(term  ^  labs- 
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This  implies  that  the  rate  of  the  termination  reaction  has  essentially  no  effect  on  the 
total  recombination  probability,  since  ')rec  ~  ‘^'labs^  and  therefore  does  not  influence 
any  of  the  particle  or  heat  fluxes. 

The  rate  of  the  abstraction  reaction  (and  thus  of  hydrogen  recombination  at  the 
diamond  surface),  on  the  other  hand,  has  strong  and  complementary  affects  on  the 
hydrogen-atom  and  recombinative  heat  fluxes  in  the  low-  and  high-pressure  regimes. 
At  low  pressures,  the  flux  of  H  atoms  is  essentially  independent  of  7rec;  thus  the  re- 
combinative  heat  flux  (which  is  just  proportional  to  7rec</'H)  increases  linearly  with  the 
recombination  probability.  At  high  pressures,  the  rate  of  replenishment  of  hydrogen 
atoms  at  the  substrate  is  limited  by  their  rate  of  diffusion  from  the  activating  surface 
through  the  H2  population.  This  (fixed)  replenishment  rate  balances  the  loss  rate 
7rec0H  of  H  atoms  at  the  diamond  surface.  Thus,  the  total  H-atom  flux  decreases 
inversely  with  the  recombination  probability;  and  the  associated  recombinative  heat 
flux  is  then  independent  of  ')rec- 

In  contrast  to  these  effects  on  <f)y{  and  qreci  the  abstraction  and  recombination 
probabilities  have  a  negligible  impact  on  the  conductive  heat  flux  qcond-  At  low  pres¬ 
sures,  the  composition  of  the  gas  is  sensitive  to  the  relative  rates  of  H-atom  recombi¬ 
nation  and  H2-moIecule  dissociation,  but  the  thermal  conductivities  of  the  two  gases 
are  very  similar.  At  high  pressures,  the  difference  between  the  conductivities  of  the 
gases  is  greater,  but  the  composition  profiles  are  only  weakly  dependent  upon  the 
recombination  probability.  Thus,  the  conductive  heat  flux  varies  very  little  with  7rec 
over  the  entire  pressure  range. 

Finally,  our  standard  conditions  assume  that  the  H2  dissociation  probability  at 
the  activating  surface  is  unity.  Generally,  the  gas  composition  and  the  particle  and 
heat  fluxes  are  essentially  independent  of  the  H2  dissociation  rate  so  long  as  it  is  much 
faster  than  the  H  recombination  rate.  When  those  two  rates  become  comparable,  the 
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H-atom  and  recombinative  heat  fluxes  decline  at  low  pressures  due  to  the  decrease 
in  the  H-atom  concentration.  These  effects  gradually  fade  at  higher  pressures,  as 
the  species  profiles  become  independent  of  the  reaction  probabilities.  The  conductive 
heat  flux  is  essentially  unaffected  by  'ydus  at  all  pressures. 

Our  numerical  model  can  be  extended  and  improved  to  make  it  more  applicable 
to  real  diamond  CVD  reactors.  For  example,  in  plasma-enhanced  CVD  (PECVD), 
molecular  hydrogen  dissociation  and  heating  of  the  gas  occurs  through  electron- 
stimulated  reactions  rather  than  interactions  with  an  activating  surface.  Using  the 
appropriate  collision  cross  sections  and  electron  densities,  the  surface  heating  and 
dissociation  can  be  replaced  by  volumetric  processes.  In  the  present  model,  carbon 
addition  to  the  diamond  surface  is  ignored.  Future  work  should  treat  the  growth 
species  as  a  perturbation  on  the  pure  H/H2  mixtures,  in  which  case  diamond  growth 
mechanisms  also  could  be  introduced  into  the  model. 

In  conclusion,  we  have  shown  that  the  DSMC  method  can  be  used  to  both  predict 
general  trends  and  obtain  quantitatively  accurate  results  for  the  heat  and  hydrogen- 
atom  fluxes  to  the  surface  in  diffusion-dominated  CVD  diamond  reactors.  While 
excellent  agreement  was  found  between  the  numerical  and  analytical  approaches  in 
the  low-  and  high-pressure  limits,  the  DSMC  technique  has  its  greatest  utility  in  the 
analytically  intractable  regime  between  these  two  regimes. 
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APPENDIX:  DERIVATION  OF  ANALYTIC  FORMS 


1.  Low-pressure  regime 

At  low  pressures,  the  counterstreaming  populations  of  particles  from  the  acti¬ 
vating  surface  and  the  substrate  are  assumed  to  be  thermalized  distributions  (i.e., 
one-sided  Maxwellians)  at  the  temperatures  Tact  ^'iid  respectively.  Denoting  the 
scalar  velocity  component  normal  to  the  wall  ux  and  the  transverse  vector  component 
V||,  the  velocity  distribution  is 

{O  (OU  'T/  ^3/2  l2kBT 

2n  (2fcBT/7rm)  /  e  V 

0, 

The  particle,  momentum,  and  energy  fluxes  associated  with 
to  be 


otherwise. 

(Al)  are  readily  shown 


<f>±  —  {^kgTI'KTny'^^  n. 

(A2) 

Mx  =  nksT, 

(A3) 

Q±  =  {SksT nksT, 

(A4) 

respectively. 

Let  and  be  the  number  densities  of  particles  of  species  s  leaving  the 


29 


activating  surface  and  the  substrate,  respectively.  For  a  nonreactive  species  such  as 
CO,  equating  the  particle  fluxes  in  the  two  directions  gives,  using  (A2), 


T, 


1/2  ^ 
act 


n 


sub 

CO 


rpl/2 

sub  ' 


(A5) 


Given  a  total  number  density  nco»  the  hot  and  cold  populations  are 


^act  _ 

^co  — 


pl/2 

sub 


-.1/2 


+  r, 


1/2 


^sub 

^co 


rpl/2 

^  act 


1/2 

sub 


nco, 


nco- 


(A6) 


The  corresponding  relations  for  the  reactive  species  H  and  H2  are  only  slightly  more 
complex.  Of  the  H  atoms  streaming  toward  the  substrate,  a  fraction  7rec/2  undergo 
abstraction  reactions,  rebounding  as  H2  molecules,  and  an  equal  fraction  participate 
in  termination  reactions,  sticking  to  the  substrate.  Thus,  the  remaining  fraction 
1  —  Tree  of  the  H  atom  flux  streaming  toward  the  substrate  must  be  balanced  by  the 
flux  leaving  that  surface,  viz., 


(1  -  7rec) 


rpll2  _  ^Sub  rpl/2 

act  —  ^  sub  ’ 


(A7) 


whence 


^act 

nn  = 


rpl/2 


sub  _  (1  ~  Irec)  T^ct 


nn, 


nn. 


(A8) 


V  ^  Krec )  act  *  J-  sub 

Similarly,  the  non-dissociating  fraction  1  —  ')diss  of  the  H2  flux  approaching  the  acti¬ 


vating  wall  balances  the  flux  leaving  that  surface, 


(A9) 


so  that 
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(AlO) 


Substituting  from  (A8)  and  introducing  the  fractional  concentrations  Xg  =  Usl'^s'^s’ 
of  species  s  yields  Eq.  (6)  of  §IV. 

The  relative  concentrations  of  hydrogen  atoms  and  molecules  are  fixed  by  noting 
that  the  mass  flux  of  H  atoms  converted  to  either  H2  molecules  or  terminated  sites  at 
the  substrate  must  balance  the  mass  flux  of  H2  molecules  dissociating  at  the  activating 
surface, 

Irecn^^'  {mnTactf''  =  •  (A12) 


Eq.  (10)  follows  upon  substitution  of  (A8)  and  (AlO)  into  this  result.  The  gas  pressure 
is  determined  by  computing  the  total  momentum  fluxes  at  either  the  substrate  or  the 
activating  wall.  From  (A3), 

P  =  +  n%%)  keTaa  +  +  nS*/  +  n^“^)  (A13) 

which  reduces  to  Eq.  (7)  upon  use  of  (A6),  (A8),  and  (AlO).  Finally,  the  energy 
flux  to  the  substrate  and  from  the  activating  surface  are  computed  from  (A4),  with 
two  modifications.  First,  the  excited  internal  degrees  of  freedom  of  the  diatomic 
molecules  must  be  taken  into  account.  At  the  temperatures  of  interest  in  diamond 
CVD,  the  rotational  modes  are  fully  excited  while  the  vibrational  modes  are  only 
weakly  populated.  This  enhances  the  heat  flux  (A4)  by  a  contribution  (f>xkBT,  and 
thus  a  multiplicative  factor  of  3/2,  for  H2  and  CO.  For  the  conductive  heat  flux  we 
thus  obtain 
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q=  +  jnJ',5’"co  ’ 

-  nS^m5''“+ jn£‘mS;'"  +  |nSgmci''"](8/ir)’'"(/CBr„6f'".  (A14) 

Substituting  from  (A6),  (A8),  and  (AlO)  and  rearranging  using  (All)  and  (A12) 
yields  Eq.  (8)  of  §IV.  Second,  the  heat  of  formation  of  H  atoms  from  H2  molecules 
must  be  taken  into  account.  The  resulting  recombinative  heat  flux  is  just  the  product 
of  the  particle  flux  the  fraction  7rec  of  recombining  H  atoms,  and  the  binding 
energy  per  atom,  ksTtindl^, 

Qrec  ~  4^H'yreckBTl)indl^t 

as  given  in  Eq.  (9)  of  §IV. 

2.  High-pressure  regime 

At  high  gas  pressures,  the  collisional  mean  free  path  of  the  hydrogen  atoms  and 
molecules  is  sufficiently  short  compared  to  the  separation  between  the  substrate  and 
the  activating  surface  that  a  hydrodynamic  description  becomes  valid.  The  contin¬ 
uum  transport  equations  that  result  are  highly  nonlinear  and  generally  require  the 
use  of  numerical  methods  to  obtain  a  complete  solution.  Such  an  approach  is  much 
too  ambitious  for  our  present  purpose.  We  will  content  ourselves  with  an  approx¬ 
imate  analysis  of  the  H/H2  system,  which  nevertheless  yields  quantitatively  useful 
predictions  as  well  as  the  correct  qualitative  trends. 

Two  features  of  the  high-pressure  solutions  can  be  explained  by  recourse  to  sim¬ 
ple  physical  arguments.  Because  the  collisions  in  a  gas  are  predominantly  binary 
interactions,  the  mean  free  path  of  any  gas  particle  -  the  distance  over  which  it  can 
transport  its  mass,  momentum,  or  energy  -  varies  inversely  with  the  number  density 
of  scattering  atoms  and  molecules.  At  the  same  time,  the  population  of  transporting 
particles  is  directly  proportional  to  the  number  density.  The  mass,  momentum,  or 
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energy  flux  carried  is  determined  by  the  product  of  these  two  factors,  which  is  there¬ 
fore  independent  of  the  number  density,  he.,  of  the  pressure.  This  accounts  for  the 
observed  saturation  of  the  H-atom  and  heat  fluxes  at  high  pressures  P .  These  fluxes 
are  driven  by  gradients  in  the  macroscopic  variables  (density,  temperature,  composi¬ 
tion,  etc.).  In  the  hydrodynamic  regime,  the  gradients  are  relatively  shallow,  so  that 
the  fluxes  are  linear  in  the  gradients  and  thus  vary  as  the  inverse  first  power  of  the 
macroscopic  scale  lengths.  This  accounts  for  the  observed  l/L  scaling  of  the  fluxes 
at  high  pressures. 

In  a  collisional  gas  mixture,  each  species  s  can  acquire  a  macroscopic  diffusion 
velocity  v,.  These  velocities  are  driven  by  the  concentration  gradients  at  a  rate 
determined  by  the  binary  diffusion  coefficients  Ds-s>  of  species  s  through  species  s'. 
They  satisfy  the  Stefan-Maxwell  equations  [29] 


dX.  . 

=  Z^— — -n{vs>  -  Vs), 


(A16) 


dx 

where  the  coefficient  ds-s>  =  nDs-s'  generally  depends  upon  the  temperature  but 
is  independent  of  the  number  density.  For  a  binary  mixture  of  H  atoms  and  H2 
molecules,  Eq.  (A16)  yields 


^  =  -i-  (XhtiXh.uh,  -  XH,nXHUH)  .  (A17) 

ox  ClH-Hz 

In  steady  state,  the  particle  fluxes  must  be  independent  of  position  and  the  net  mass 
flux  must  be  zero,  whence 


nX\{Vi{  =  —2nX}i2Vii2  —  —  constant. 


(A18) 


where  ip  is  the  magnitude  of  the  net  H  atom  flux  toward  the  substrate  (at  x  =  0). 
Substituting  (A18)  into  (A17)  and  using  Xn^  =  1  ~  gives 


dXn 

dx 


tp 

2dH-H2 


(Ah -2), 


(A19) 
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whose  solution  is 


Xh  =  2  -  ylexp(-xi/)/2dH-H2)-  (A20) 

We  have  neglected  the  temperature  dependence  of  dH-H2  here,  and  A  is  a  constant 
to  be  determined. 

It  remains  to  apply  boundary  conditions  to  the  species  concentrations  Xs-  At  suf¬ 
ficiently  high  pressures,  the  populations  adjacent  to  the  substrate  and  the  activation 
surface  will  be  completely  depleted  of  H  atoms  (which  recombine)  and  H2  molecules 
(which  dissociate),  respectively.  In  the  limit,  then,  we  have 

=  A:h(0)  ^  0  =;>  A  =  2,  (A21) 

X^^  =  Xn{L)^l  xP  =  \nidn.nJL.  (A22) 

It  is  crucial  to  note  that  V’  is  the  net  flux  of  hydrogen  atoms  from  the  activating 
surface  to  the  substrate.  The  total  flux  of  H  atoms  striking  the  substrate,  (/>h,  is  just 
the  particle  flux  associated  with  an  isotropic  Maxwellian  distribution,  viz., 

<f>u  =  {kBTsub/2irrn„f^  =  X^^^P  (27rmHfcBT,„6)"'/'* .  (A23) 

The  net  H-atom  flux  0  must  balance  that  fraction  7rec  of  the  total  flux  (^h  which 
recombines  at  the  substrate, 

4’  =  'yrec<l>n-  (A24) 

Combining  (A24)  and  (A22),  and  averaging  the  diffusion  coefficient  at  the  substrate 
and  activating  temperatures  to  approximately  account  for  its  temperature  depen¬ 
dence,  we  obtain 

»  ln2  +  dJi'in,)  h.„L  (A25) 

for  the  total  atomic  hydrogen  flux  to  the  substrate. 

Results  (A23)  and  (A25)  can  be  combined  to  yield 
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«  {2xmnkBT,ubf  (A26) 

Tree  ^ 

Thus,  <x  P~^  -+  0  as  the  pressure  increases,  as  asserted  above.  The  net  Hj  flux 
to  the  activating  surface,  i/’/2,  is  related  to  the  total  flux  (/1H2  t>y 

V>/2  =  7dis3<^H2  5  (^27) 

whence 

XtS  «  (A28) 

Consequently,  ->  0  at  high  pressures,  as  well. 

The  total  energy  flux  is  composed  of  two  terms:  the  familiar  conductive  compo¬ 
nent  given  by  Fourier’s  law,  i.e.,  the  product  of  the  temperature  gradient  and  the 
thermal  conductivity  of  the  gas;  and  a  contribution  due  to  the  flow  of  enthalpy  at  the 
species’  diflFusion  velocities.  Thus, 

=  +  (A29) 

where  the  specific  enthalpies  are  h\{  =  ksTbindj"^  +  ^ksT  12  and  =  IksT (2  for 
atomic  and  molecular  hydrogen,  respectively.  The  recombinative  heat  flux  computed 
in  the  DSMC  model  is  just  the  contribution  due  to  the  heat  of  formation  of  H  from 
H2,  whence 

9rec  =  <t>H'llreckBTbindl^-  (A30) 

The  remaining  terms  in  (A29)  constitute  the  conductive  heat  flux  as  calculated  by 
the  DSMC  model.  For  a  gas  mixture  such  as  H/H2,  the  thermal  conductivity  k  is  a 
summation  of  the  conductivities  of  the  individual  constituents.  Each  such  contribu¬ 
tion  is  weighted  by  that  species’  concentration  and  reduced  by  a  factor  accounting  for 
binary  collisions  with  other  species  present  in  the  mixture  [29].  Also,  the  individual 
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conductivities  vary  with  temperature.  We  avoid  these  complexities  by  using  the  mean 
of  the  conductivities  at  the  substrate  and  activating  temperatures,  and  the  average 
temperature  gradient  across  the  system,  to  estimate  the  heat  flux.  As  noted  above, 
at  high  pressures  the  composition  tends  toward  100%  H  atoms  at  the  activating  sur¬ 
face  and  100%  H2  molecules  at  the  substrate.  Then  the  conductive  heat  flux  can  be 
approximated  as 

qcond  «  (kh*  +  -  Tsub)  /2T  +  3<^H7rec^S  {Tact  +  Tsub)  /8,  (A31) 

as  given  in  Eq.  (16)  in  §4.  Here,  we  have  used  the  average  of  the  activating  and 
substrate  temperatures  to  estimate  the  enthalpic  contribution  ^H7rec3A:jgr/4  to  the 
conductive  flux. 
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TABLES 


TABLE  1.  Parameters  for  the  variable  soft-sphere  collision  model:  d  is  the  molec¬ 
ular  diameter,  the  reference  temperature,  uj  the  viscosity  index,  a  the  scattering 
parameter,  m  the  molecular  mass,  and  Z’’®*  the  rotational  relaxation  number.  Values 
were  taken  from  Bird  [3]  where  available.  The  values  of  u;  and  a  for  atomic  hydrogen 
were  set  equal  to  those  for  helium.  The  diameter  for  atomic  hydrogen  and  the  rota¬ 
tional  relaxation  numbers  for  molecular  hydrogen  and  carbon  monoxide  were  taken 
from  the  Sandia  transport  data  base  [19].  The  parameters  Ci  and  C2  for  the  temper¬ 
ature  dependence  of  the  vibrational  relaxation  number  (Z„  =  {Ci  fT'^)  exp{C2T~^^^)) 
were  set  to  9.1  and  220.0,  respectively,  for  all  collisions  involving  H2. 


species 

d(k) 

T-/  (K) 

U 

a 

m  (amu) 

Z"®‘(298) 

H2 

2.88 

273 

0.67 

1.35 

2 

280.0 

CO 

4.12 

273 

0.73 

1.49 

28 

1.8 

H 

2.05 

273 

0.66 

1.26 

1 

— 

TABLE  2.  Comparison  between  free-streaming  analytic  (indicated  by  superscript 
FS)  and  low-pressure  DSMC  results  for  H/H2  mixtures.  The  conductive  (qcond)  and 
recombinative  {qrec)  heat  fluxes  are  in  W  cm”^.  L  is  the  separation  between  the 
substrate  and  activating  surface,  P  is  the  pressure,  and  Ah  is  the  atomic  hydrogen 
mole  fraction.  was  evaluated  by  averaging  Ah  at  the  bounding  surfaces.  In 

all  cases,  the  spatial  variation  in  Ah  was  less  than  2%. 
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L  (cm) 

P  (Torr) 

DSMC 

Hcond 

& 

DSMC 

Hrec 

FS 

Hrec 

0.25 

0.043 

1.63 

1.65 

0.89 

0.89 

0.965 

0.967 

0.5 

0.045 

1.65 

1.73 

0.95 

0.94 

0.965 

0.967 

1.0 

0.050 

1.74 

1.93 

1.09 

1.04 

0.960 

0.967 

TABLE  3.  Comparison  between  analytic  hydrodynamic  (indicated  by  superscript 
HD)  and  high-pressure  DSMC  results  for  H/H2  mixtures.  The  conductive  {qcond)  and 
recombinative  (^rec)  heat  fluxes  are  in  W  cm“^,  and  the  atomic  hydrogen  fluxes  to 
the  substrate  are  in  10^®  atoms  cm~^  s~^.  L  is  the  separation  between  the  substrate 
and  activating  surface  and  P  is  the  pressure. 


L  (cm) 

P  (Torr) 

DSMC 

Hcond 

(lc<^d 

DSMC 

Hrec 

HD 

Hrec 

<j>DSMC 

0.25 

26.2 

90.1 

128. 

279. 

400. 

1052. 

1520. 

0.5 

25.0 

48.3 

56. 

175. 

200. 

660. 

760. 

1.0 

24.2 

25.9 

28. 

101. 

100. 

381. 

380. 
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Tsub  ^act 


H-atom  recombination  H2  dissociation 


FIG.  1.  Computational  domain  used  for  the  DSMC  calculations. 
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FIG.  2.  (a) 
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FIG.  2.  Spatially  resolved  (a)  temperature  and  (b)  species  profiles  for  H/H2 
mixtures.  The  diamond  substrate  at  1200  K  is  located  at  a:  =  0,  and  the  activating 
surface,  located  at  x  =  1  cm,  is  maintained  at  3000  K.  Surface  reaction  probabilities 
are  those  of  Brenner  et  al.  and  100%  hydrogen  dissociation  at  the  activating  surface 
is  assumed. 


40 


FIG.  2.  (b) 


FIG.  2.  Spatially  resolved  (a)  temperature  and  (b)  species  profiles  for  H/H2 
mixtures.  The  diamond  substrate  at  1200  K  is  located  at  a;  =  0,  and  the  activating 
surface,  located  at  x  =  1  cm,  is  maintained  at  3000  K.  Surface  reaction  probabilities 
are  those  of  Brenner  et  al.  and  100%  hydrogen  dissociation  at  the  activating  surface 


is  assumed. 
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FIG.  3.  (a) 


-o-  L  =  0.25  cm 
-o-  L  =  0.5  cm 
-A-  L  =  1 .0  cm 


E  200 


CT  150 


1400 


1200 


1000-e- 

I 


o 

800  s 
0) 
o 
3 

600  w 
o 
3 


gi-s — 

. *A’*' 


- - - 


0  5  10  15  20  25 

P  (Torr) 

FIG.  3.  Pressure  dependence  of  qcond  (dotted  lines),  grec  (dashed  lines),  and  qtot 
(solid  lines)  for  H/H2  mixture  at  various  values  of  the  substrate-heat  source  separa¬ 
tion.  The  hydrogen  atom  flux,  <^h  (dashed  lines),  is  exactly  proportional  to  qrec  and 
is  plotted  against  the  right  axis.  Figure  3b  is  a  blow  up  of  the  low-pressure  region  of 

g.  3a. 


FIG.  3.  (b) 
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FIG.  3.  Pressure  dependence  of  qcond  (dotted  lines),  qrec  (dashed  lines),  and  qtot 
(solid  lines)  for  H/H2  mixture  at  various  values  of  the  substrate-heat  source  separa¬ 
tion.  The  hydrogen  atom  flux,  (dashed  lines),  is  exactly  proportional  to  qrec  and 
is  plotted  against  the  right  axis.  Figure  3b  is  a  blow  up  of  the  low-pressure  region  of 

Fig.  3a. 
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(1)^  (10^®  atoms/cm^-s) 
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FIG.  4.  (a) 


FIG.  4.  Spatially  resolved  (a)  temperature  and  (b)  species  profiles  for  H/H2/CO 
mixtures.  The  diamond  substrate  at  1200  K  is  located  at  x  —  0,  the  activating 
surface,  located  at  x  =  1  cm,  is  maintained  at  3000  K.  Surface  reaction  probabilities 
are  those  of  Brenner  et  al.  and  100%  hydrogen  dissociation  at  the  activating  surface 
is  assumed. 
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FIG.  4.  (b) 


mixtures.  The  diamond  substrate  at  1200  K  is  located  at  a;  =  0;  the  activating 
surface,  located  at  a;  =  1  cm,  is  maintained  at  3000  K.  Surface  reaction  probabilities 
are  those  of  Brenner  et  al.  and  100%  hydrogen  dissociation  at  the  activating  surface 
is  assumed. 
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FIG.  5.  Comparison  of  the  pressure  dependence  of  qcond  (dotted  lines),  qrec 
(dashed  lines),  and  qtot  (solid  lines)  for  H/H2  and  H/H2/CO  mixtures.  L  is  1  cm.  The 
hydrogen  atom  flux,  <^h  (dashed  lines),  is  exactly  proportional  to  qrec  and  is  plotted 
against  the  right  axis.  Figure  5b  is  a  blow  up  of  the  low-pressure  region  of  Fig.  5a. 
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FIG.  5.  Comparison  of  the  pressure  dependence  of  qcond  (dotted  lines),  qrec 
(dashed  lines),  and  qtot  (solid  lines)  for  H/H2  and  H/H2/CO  mixtures.  L  is  1  cm.  The 
hydrogen  atom  flux,  (^h  (dashed  lines),  is  exactly  proportional  to  qrec  and  is  plotted 
against  the  right  axis.  Figure  5b  is  a  blow  up  of  the  low-pressure  region  of  Fig.  5a. 
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(10^®  atoms/cm^-s) 


FIG.  6.  (a) 


FIG.  6.  Pressure  dependence  of  qcond  (dotted  lines),  qrec  (dashed  lines),  and  qtot 
(solid  lines)  for  H/H2  mixtures  for  various  values  of  the  temperature  at  the  activating 
surface.  L  is  1  cm.  The  hydrogen  atom  flux,  (dashed  lines),  is  exactly  proportional 
to  qrec  and  is  plotted  against  the  right  axis.  Figure  6b  is  a  blow  up  of  the  low-pressure 
region  of  Fig.  6a. 
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FIG.  7.  Pressure  dependence  of  qcond  (dotted  lines),  g^ec  (dashed  lines),  and  qtoi 
(solid  lines)  for  H/H2  mixtures  at  various  values  of  the  hydrogen  dissociation  proba¬ 
bility.  L  is  1  cm.  The  hydrogen  atom  flux,  <^h  (dashed  lines),  is  exactly  proportional 
to  qrec  and  is  plotted  against  the  right  axis.  Figure  7b  is  a  blow  up  of  the  low-pressure 


region  of  Fig.  7a. 
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FIG.  7.  (b) 


P  (Torr) 

FIG.  7.  Pressure  dependence  of  gcond  (dotted  lines),  grec  (dashed  lines),  and  qtot 
(solid  lines)  for  H/H2  mixtures  at  various  values  of  the  hydrogen  dissociation  proba¬ 
bility.  L  is  1  cm.  The  hydrogen  atom  flux,  (f>}i  (dcished  lines),  is  exactly  proportional 
to  qrec  and  is  plotted  against  the  right  axis.  Figure  7b  is  a  blow  up  of  the  low-pressure 
region  of  Fig.  7a. 
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FIG.  8.  Spatially  resolved  species  profiles  for  H/H2  mixtures  with  hydrogen  dis¬ 
sociation  probabilities  of  100%  and  10%  at  the  activating  surface.  The  diamond 
substrate  at  1200  K  is  located  at  x  =  0;  the  activating  surface,  located  at  x  =  1  cm, 
is  maintained  at  3000  K. 
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1  cm.  The  hydrogen  atom  flux,  (dashed  lines),  is  exactly  proportional  to  3'>id 
is  plotted  against  the  right  axis.  Figure  9b  is  a  blow  up  of  the  low-pressure  region  of 


Fig.  9a. 
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FIG.  9.  Pressure  dependence  of  qcond  (dotted  lines),  ^rec  (dashed  lines),  and  qtot 
(solid  lines)  for  H/H2  mixtures  for  three  sets  of  the  surface  reaction  probabilities.  L  is 
1  cm.  The  hydrogen  atom  flux,  (dashed  lines),  is  exactly  proportional  to  qrec  and 
is  plotted  against  the  right  axis.  Figure  9b  is  a  blow  up  of  the  low-pressure  region  of 


Fig.  9a. 


